Van der Waals Interactions Between Thin Metallic Wires and Layers 
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Quantum Monte Carlo (QMC) methods have been used to obtain accurate binding-energy data 
for pairs of paraUel, thin, metaUic wires and layers modelled by ID and 2D homogeneous electron 
gases. We compare our QMC binding energies with results obtained within the random phase 
approximation, finding significant quantitative differences and disagreement over the asymptotic 
behavior for bilayers at low densities. We have calculated pair-correlation functions for metallic 
biwire and bilayer systems. Our QMC data could be used to investigate van der Waals energy 
functionals. 
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One-dimensional conductors such as carbon nanotubes 
are essential components of many proposed nanoteclmo- 
logical devices, and are currently the subject of numerous 
experimental and theoretical studies. It has recently been 
demonstratedi that the van der Waals (vdW) interaction 
between pairs of distant, parallel, thin conducting wires 
assumed in many current models of metallic carbon nan- 
otubes is qualitatively wrong. In this letter we provide 
the first accurate binding-energy data for pairs of thin 
metallic wires and layers, which can be used as a bench- 
mark for subsequent theoretical studies or to parametrize 
model interactions between ID and 2D conductors. 

Thin, electrically neutral wires are attracted to one 
another by vdW forces. The standard method for cal- 
culating the vdW interaction between objects is to as- 
sume there are pairwise interactions between volume ele- 
ments with an attractive tail of the form Upp{r) oc —r~^, 
which is appropriate for the vdW interaction between 
molecules. Summing these interactions for a pair of ID 
parallel wires separated by a distance z gives a vdW 
binding energy of U{z) oc —z^^. Such pairwise vdW 
models have been used in studies of single-walled carbon 
nanotubes.— However, a recent investigatiort^ of the in- 
teraction between pairs of thin, metallic wires modelled 
by ID homogeneous electron gases (HEGs) within the 
random phase approximation (RPA) found that the bind- 
ing energy falls off (approximately) as^ 
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where b is the wire radius and 2rg is the length of the wire 
section containing one electron on average. The pairwise 
vdW model is clearly appropriate for an insulator or for 
a metallic wire whose radius is greater than the screening 
length, but is inappropriate for a thin conductor such as 
a single-walled carbon nanotubcfi 

Likewise, the binding energy per particle of a pair of 
thin, parallel metallic layers can be shown to decay as 
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within the RPA,— compared with U{z) oc —z ^ within 
the pairwise vdW theory, where z is the layer separation. 



(In a 2D HEG is the radius of the circle that contains 
one electron on average.) At very large separations the 
vdW attraction is dominated by the Casimir effect in 
which the zero-point energy of photon modes between the 
metallic layers gives rise to an attractive force. However 
we restrict our attention to the range of separations in 
which vdW effects are dominant. 

Within the RPA, the binding energy may be calculated 
as the change in the zero-point energy of plasmon modes 
as a function of separationii However, the RPA is poor 
in low-dimensional systems, and ceases to be valid at 
low densities, where correlation effects become dominant. 
We have therefore performed quantum Monte Carlo^i^ 
(QMC) calculations of the binding energies of pairs of 
thin, metallic wires and layers modelled by ID and 2D 
HEGs with neutralizing backgrounds. In particular we 
have used the variational and diffusion quantum Monte 
Carlo (VMC and DMC) methods as implemented in the 
CASINO code<^ DMC is the most accurate method avail- 
able for studying quantum many-body systems such as 
electron gases. We have also calculated pair-correlation 
functions (PCFs), enabling us to examine the correlation 
hole responsible for the vdW attraction between pairs of 
wires and layers. 

In our QMC calculations we use the full Coulomb 
potential, so that for ID HEGs the many-electron 
wave function must go to zero at both parallel- and 
antiparallel-spin coalescence points for electrons in the 
same wire. The nodal surfaces for paramagnetic and 
ferromagnetic ID HEGs are therefore the same, and so 
the fixed-node DMC energy — which is equal to the exact 
ground-state energy because the nodal surface is exact — 
is independent of the spin polarization. This conclusion 
does not violate the Lieb-Mattis theoremiS because the 
ID Coulomb interaction is pathological in the formal 
sense of Lieb and MattiSf^i For convenience we choose 
to work with ferromagnetic ID HEGs, because a Slater 
determinant wave function produces the correct nodal 
surface in this case. Our QMC studies of ID HEGs will 
be published elsewhereii^ 

Fermionic symmetry is imposed in 2D via the fixed- 
node approximation,— in which the nodal surface is con- 
strained to equal that of a trial wave function. We ex- 
pect a very high degree of cancellation of fixed-node er- 
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rors when the binding energy is calculated. It has already 
been shown that fixed-node DMC is able to describe vdW 
forces between heliunJ^ and neon atomsfi^ 

We have verified that time-step and population-control 
biases in our DMC energies are negligible by repeating 
some of the calculations using different time steps and 
populations. Finite-size errors arc a more serious prob- 
lem, although most of the bias cancels out when the 
energy difference is taken to obtain the binding energy. 
Twist averagingi^ or the addition of finite-size corrections 
are unlikely to reduce the bias in the binding energy. We 
expect our binding-energy results to be valid so long as 
the wire or layer separation is small compared with the 
length of the ID or 2D simulation cell; for larger sepa- 
rations the system resembles a pair of insulators and the 
binding energy is expected to fall off more steeply (in 
accordance with the pairwise vdW model) . 

We have used Slater- Jastrow-backfiow trial wave 
functions For our 2D bilayer calculations the Slater 
part of the wave function consists of a product of four 
determinants of plane- wave orbitals for spin-up and spin- 
down electrons in each of the two layers. For our ID bi- 
wire calculations the Slater part consists of a product of 
two determinants of plane-wave orbitals for the electrons 
in each wire (recall that each wire is ferromagnetic in our 
calculations). Slater determinants for a ID HEG are of 
Vandermonde form, and could therefore be rewritten as 
polynomials and evaluated in a time that scales linearly 
with system size; however other parts of the QMC al- 
gorithm such as the evaluation of the two-body Jastrow 
terms and backflow functions take up a significant frac- 
tion of the computer time, so for convenience we have 
continued to employ the usual determinant-evaluating 
and updating machinery of QMC calculations^iS 

Our Jastrow factors consist of polynomial and 
plane-wave two-body termsi^ satisfying the Kato cusp 
conditionsi^ In spite of the fact that the nodal surface is 
exact in ID, two-body backflow correlationsii were found 
to make a very significant improvement to the wave func- 
tion, as can be seen in Table HI Backflow functions were 
used in all of our ID calculations and in our 2D calcu- 
lations at = 1 a.u. Free parameters in the trial wave 
function were optimized by minimizing the unreweighted 
variance of the local energyj ^^i^^ 

Each wire or layer is accompanied by a neutralizing 
background. When studying biwircs (bilayers) wc need 
to add the interaction between the electrons in each wire 
(layer) and the background of the opposite wire (layer), 
plus the interaction of the two backgrounds. This contri- 
bution to the energy is i?cap = log(z)/(2rs) for a biwirc 
and ii'cap = z/r1 for a bilayer. 

The binding energies of pairs of ID and 2D HEGs arc 
shown in Figs. [1] and (2] The approximate RPA binding 
energy shown in Fig. [1] was obtained using b = rs/10 
in Eq. ([T]). The wire radius h is therefore small com- 
pared with the other length scales in the system. We 
have included DMC results at several different system 
sizes in Figs. [T] and [H In ID, half the length of the sim- 



Method 


Energy (a.u. / el.) 


Var. (a.u.) 


%age corr. en. 


HF 


-0.215943040112 




0% 


SJ-VMC 


-0.2319668(4) 


0.0000360(2) 


99.974(3)% 


SJB-VMC 


-0.2319710(3) 


0.0000046(1) 


100.000(3)% 


SJB-DMC 


-0.2319709(3) 




100% 



TABLE I; Energy, energy variance, and fraction of corre- 
lation energy retrieved using different levels of theory and 
wave function for a 15-electron ID ferromagnetic HEG at 
rs=15 a.u. "HF" stands for Hartree-Fock theory, "SJ" de- 
notes a Slater- Jastrow trial wave function and "SJB" a Slater- 
Jastrow-backflow trial wave function. 

ulation cell is L^^.n = Nr^^ where N is the number of 
electrons per wire, and in 2D the size of the simulation 
cell is Lr^^N ~ VNts- The binding energy falls off more 
steeply once z becomes a significant fraction of L,.^ jv, as 
expected. Clearly the binding energies enter the asymp- 
totic regime when z ^ r^. We have therefore fitted the 
RPA asymptotic binding-energy forms to our QMC data 
in the range ^ z ^ L^^^n- We believe the errors in 
the fitted exponents are about 0.1-0.2. 

The fits to the DMC biwire binding-energy data shown 
in Fig. [T] are 

Ui{z) = -0.0815Z-2-28 [log(27000z)]"^/^ (3) 
Uaiz) = -0.0225^-1-^* [log(1.95z)]"^/' (4) 
t/io(z) = -0.0967z"2■l^[log(0.4922)]~^/^ (5) 

where Ur^ {z) is the binding energy at density parameter 
Tg. The DMC binding-energy data are clearly in much 
better agreement with the RPA [Eq. ^] than with the 
pairwise vdW theory [U{z) cx z~^]. It is not meaningful 
to compare the prefactors because of the arbitrariness of 
our choice of the wire radius b in the RPA theory. 

The fits to the DMC bilayer binding-energy data shown 
in Fig. [5] arc 

Ui{z) = -0.00637z-2-5*^ (6) 
Uaiz) = -OmSSz-^-^'^ (7) 
Uioiz) = -0.882z-3-i6^ (8) 

where Ur^ (z) is the binding energy at density parame- 
ter rg. At high densities (ts = 1 and 3 a.u.) our results 
clearly show the — behavior predicted by Eq. ([2]). 
At low density (r^ = 10 a.u.) the binding energy falls 
off more steeply than predicted by the RPA, although 
the asymptotic behavior is clearly better described by 
the RPA than the pairwise vdW theory [U{z) cx — z~^]. 
DMC and the RPA give similar prefactors for the asymp- 
totic binding energy at = 3 a.u., but the DMC prefac- 
tor is somewhat lower at Ts = 1 a.u. 

PCFs were accumulated by binning the interparticle 
distances in the electron configurations generated by the 
VMC and DMC algorithms. The error in the VMC and 
DMC PCFs (?vMC a-iid gnuc is first order in the error in 
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FIG. 1: (Color online) Binding energy per particle of a ID 
HEG biwire as a function of wire separation for rs = 1 a.u. 
(top panel), rs = 3 a.u. (middle panel), and = 10 a.u. 
(bottom panel). The DMC time steps were 0.04, 0.2, and 
2.5 a.u. at = 1, 3, and 10 a.u., and the target configuration 
population was 2048 in each case. 



FIG. 2: (Color online) Binding energy per particle of a 2D 
HEG bilayer as a function of layer separation for rs = 1 a.u. 
(top panel), rs = 3 a.u. (middle panel), and rs = 10 a.u. 
(bottom panel). The DMC time steps were 0.007, 0.05, and 
0.5 a.u. at rs = 1, 3, and 10 a.u., and the target configuration 
populations were 1024, 320, and 1024. 
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FIG. 3: (Color online) VMC PCF for a ID HEG biwire at 
Vs — 3 a.u. and wire separation z — 3 a.u. The inset shows 
the interwire correlation hole in greater detail. It was verified 
for smaller system sizes that the DMC and VMC PCFs were 
in excellent agreement (as expected due to the accuracy of 
the trial wave function illustrated in Table 
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the trial wave function, but the error in the extrapolated 
PCF (7oxt = 2(7dmc -~ 5VMC is second order in the error 
in the wave function.— PCFs for biwires and bilayers are 
sliown in Figs. [3] and 21 respectively. The correlation 
holes between the electrons in opposite wires or planes are 
exceedingly shallow, although they extend over a distance 
roughly equal to the separation. It can be seen that the 
interwire PCF has the long-ranged oscillatory behavior 
exhibited by ID HEG PCFs. 



In conclusion we have used QMC to obtain the first 
accurate binding-energy data for pairs of thin, parallel, 
metallic wires and layers. Our results are in broad agree- 
ment with recent RPA calculations of the binding energy 
and complete disagreement with the standard pairwise 
vdW model. However there are significant differences 
between the DMC and RPA results for bilayers: at high 
densities the asymptotic behavior of the binding energy 
as a function of separation is the same but the prefactor is 
different, and at low densities the DMC binding energy 
falls off more rapidly, implying that correlation effects 
neglected in the RPA are important in this regime. Our 
data can serve as a benchmark for future theoretical stud- 
ies of the binding energies of ID and 2D HEGs, and can 
also be used to parametrize model interactions between 
thin conductors. Together with our results for biwire and 
bilayer PCFs, our data could be used to investigate en- 
ergy functionals that incorporate vdW effects for use in 
density- functional calculations . 
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FIG. 4: (Color online) Extrapolated PCF for a 2D HEG bi- 
layer at rs = 1 a.u. and layer separation z = 2 a.u. The inset 
shows the interwire correlation hole in greater detail. The 
VMC and DMC PCFs are in excellent agreement, implying 
that the extrapolated PCF is reliable. 
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